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1. Introduction 



Gravitational microlensing (in the local group) refers to the temporal brightening of a back- 
ground star due to intervening objects. Einstein (1936) first examined (micro)lensing by a single 
star, and concluded that "there is no great chance of observing this phenomenon." Some important 
works were performed in intervening years by [46] and [34], but the research topic was revitalised 
by Paczyriski (1986) who proposed it as a method to detect compact dark matter objects in the 
Galactic halo. 

The original goal is now out of favour, since we know with high precision that most of the 
dark matter must be non-baryonic, e.g. from observations of microwave background radiation and 
nucleosynthesis (at the time of his paper, this was, however, unclear). Nevertheless, gravitational 
microlensing has turned into a powerful technique with diverse applications in astrophysics, includ- 
ing the study of the structure of the Milky Way, stellar atmospheres and the detection of extrasolar 
planets and stellar-mass black hole candidates. The field has made enormous progress in the last 



two decades. There have been a number of reviews on this topic (e.g. [g9|, |3g, |20|, |54p), the most 



recent highlight was given in [26]. This article gives an introduction to microlensing, aimed at a 
level for a starting PhD student. Together with other talks in the workshop and proceedings ' , one 
can gain a thorough feeling about the state-of-the-art research in this field (as of 2008). 

The reference list given here is seriously incomplete (and biased). For more complete refer- 
ences and information about ongoing microlensing surveys, see the review papers mentioned above 
and the web site: http : / /miens . net/ (built by Szymon Kozlowski, Subo Dong and Lukasz 
Wyrzykowski). 



2. What is gravitational microlensing? 

The light from a background source is deflected, distorted and (de)magnified by intervening 
objects along the line of sight. If the lens, source and observer are sufficiently well aligned, then 
strong gravitational lensing can occur. Depending on the lensing object, strong gravitational lensing 
can be divided into three areas: microlensing by stars, multiple-images by galaxies, and giant arcs 
and large-separation lenses by clusters of galaxies. For microlensing, the lensing object is a stellar- 
mass compact object (e.g. normal stars, brown dwarfs or stellar remnants [white dwarfs, neutron 
stars and black holes]); the image splitting in this case is usually too small (of the order of milli- 
arcsecond in the local group) to be resolved by ground-based telescopes, thus we can only observe 
the magnification change as a function of time. 

The left panel in Fig. |l] illustrates the microlensing geometry. A stellar-mass lens moves 
across the line of sight towards a background star. As the lens moves closer to the line of sight, its 
gravitational focusing increases, and the background star becomes brighter. As the source moves 
away, the star falls back to its baseline brightness. If the motions of the lens, the observer and 
the source can be approximately taken as linear, then the light curve is symmetric. Since the 
lensing probability for microlensing in the local group is of the order of 10^^ (see section |5|), the 
microlensing variability usually should not repeat. Since photons of different wavelengths follow 
the same propagation path (geodesies), the light curve (for a point source) should not depend on 

'available at http : //pos . sissa . it /cgi-bin/reader / conf . cgi?conf id=54 
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(t-to)/t, 



Figure 1: The left panel shows a side-on view of the geometry of microlensing where a lens moves across 
the line of sight towards a background source. The right panel shows two light curves corresponding to two 
dimensionless impact parameters, mq = 0.1 and 0.3. The time on the horizontal axis is centred on the peak 
time to and is normalised to the Einstein radius crossing time fe- The lower the value of mq, the higher the 
peak magnification. For the definitions of uq and fE see section 4. 1 . 



the colour. The characteristic symmetric shape, non-repeatability, and achromaticity can be used 
as criteria to separate microlensing from other types of variable stars (exceptions to these rules will 
be discussed in section H~2|). 



3. Lens equation, image positions and magnifications 

To derive the characteristic light curve shape shown in the right panel of Fig. ||, we must look 
closely at the lens equation, and the resulting image positions and magnifications for a point source. 

3.1 Lens equation 

The lens equation is straightforward to derive. Figure ^illustrates a side-on view of the lensing 
configuration. Simple geometry yields 

T]+Dds5: = |-^, (3.1) 

where Dj, Dg and Dds are the distance to the lens (deflector), distance to the source and distance 
between the lens (deflector) and the source, fj is the source position (distance perpendicular to the 
line connecting the observer and the lens), ^ is the image position, and a is the deflection angle. For 
gravitational microlensing in the local group, D^s = D^— D^? Mathematically, the lens equation 
provides a mapping between the source plane to the lens plane. The mapping is not necessarily 
one-to-one. 

^For cosmological microlensing in an expanding universe, the distances should be taken as angular diameter dis- 
tances, and in general ^ D^—D^. See the review by Wambsganss in these proceedings on cosmological microlens- 
ing. 
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Figure 2: Illustxation of various distances and angles in the lens equation (eqs. 3.1 and 3.2). 



Dividing both sides of eq. (|3JJ) by Dj, we obtain the lens equation in angles 

^ + a = d, 



(3.2) 



where j8 = T] /Dg, = f /^d, and a = a x D^^/D^ is the scaled deflection angle. These angles are 
illustrated in Fig. ^ 

For an axis-symmetric mass distribution, due to symmetry, the source, observer and image 
positions must lie in the same plane, and so we can drop the vector sign, and obtain a scalar lens 
equation: 

j8 + a = 0. 



3.2 Image positions for a point lens 

For a point lens at the origin, the deflection angle is given by 

^ 4GM 1 - 



a 



C2 ^2 



and the value of the scaled deflection angle is 



a = — a 



where we have defined the angular Einstein radius as 



Oe 



Da 



-1/2 



M 



Y Dd/Ds \SkpcJ \0.3Mq 



1/2 



The lens equation for a point lens in angles is therefore 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 
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We can further simply by normalising all the angles by Oe, rs = [5 / Be, r = 6 / Oe, the above equation 

becomes^ ^ 

r, + - = r. (3.8) 
r 

For the special case when the lens, source and observer are perfectly aligned (rg = 0), due to axis- 
symmetry along the line of sight, the images form a ring ("Einstein" ring) with its angular size 



given by eq. (3.6). 



For any other source position / 0, there are always two images, their positions are given by 



(3.9) 

The '+' image is outside the Einstein radius (r+ > 1) on the same side of the source, while the '— ' 
image is on the opposite side and inside the Einstein radius (r_ < and < 1). The angular 
separation between the two images is 

A0 = 0E(r+-r_) = eEV4 + rs2. (3.10) 

The image separation is of the same order of the angular Einstein diameter when rs < 1 , and thus 
will be in general too small to be observable given the typical seeing from the ground {^^ one 
arcsecond); we can only observe lensing effects through magnification. One exception may be the 
VLT interferometer (VLTI) which can potentially resolve the two images. This may be important 
for discovering stellar-mass black holes since they have larger image separations due to their larger 
masses than typical lenses with mass ~ 03Mq (@, p3|]). 

The physical size of the Einstein radius in the lens plane is given by 



1/2 / X 1/2 



rn = D^Oe = . /l™^ . 2.2 AU ,4 ' ( ^\ ' . (3.11) 



So the size of the Einstein ring is roughly the scale of the solar system, which is a coincidence that 
helps the discovery of extrasolar planets around lenses. 

3.3 Image magnifications 

Since gravitational lensing conserves surface brightness, the magnification of an image is sim- 
ply given by the ratio of the image area and source area. For a very small source, we can consider 
a thin source annulus with angle A0 (see Fig. ^. For a point lens, this thin annulus will be mapped 
into two annuli, one inside the Einstein ring and one outside. 

The area of the source annulus is given by the product of the radial width and the tangential 
length dr^ x rsA(/). Similarly, each image area is dr x rA0, and the magnification is given by^ 

dr X rA^ r dr 



For the two images given in eq. ( [3.12 ), one finds 



ir, + ,/7JT4f ^ ^ ^ _ (^^-7^2+4)2^ ^^^^^ 



rs is not to be confused with the size of the star, which we denote as r^,. 
^This is a special case of the determinant of the Jacobian in the lens mapping, see section 4.3 
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The magnification of the '+' image is positive, while the '— ' image is negative. The former image 
is said to have positive parity while the latter negative^. The total magnification is given by 

At = |At+| + |At-l = ^V4=T' ^3-14) 
and the difference is identical to unity 

liU+|-|M-l = l- (3.15) 
We make some remarks about the total magnification and image separations: 

1 . When = 1 , ^ = 3/ ^/5 « 1 .342, corresponding to about 0.3 19 magnitude. Such a difference 
is easily observable^, and so the area occupied by the Einstein ring is usually taken as the 
lensing "cross-section." 

2. When —>■ oo, |jU^/jU_| — > rg'*, /i — > 1 +2rs^'^. The angular image separation is given by 
Ad = {r,+2r,-^)dE. 

3. High magnification events occur when rg — > 0. The asymptotic behaviours are ^ r,"^ (1 + 
3rs^/8), Ad (2 + rs^/4)0E> and dr/dr^ —>■ 1/2. The last relation implies that, at high 
magnification, the image is compressed by a factor of 2 in the radial direction (see Fig. 




Figure 3: Images of a thin annulus from rs to rs + drg by a point lens on the plane of the sky. The dashed 
line is the Einstein ring. A(j) is the angle subtended by the thin annulus. 



4. Light curve and microlensing degeneracy 

Given a source trajectory, we can easily describe the standard hght curve with a few parameters 
which suffers from the microlensing degeneracy. 

^Let us imagine two arrows for the thin annulus (see Fig. one in the radial direction and one in the tangential 
direction respectively. For the negative parity image, the corresponding tangential arrow for the image is reversed with 
respect to that in the source, while in the radial direction the arrow directions remain the same for the source and image. 
For the positive parity image, the directions of the arrows are the same for the image and the source. 

^For bright stars, the accuracy of photometry can reach a few milli-magnitudes. 
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4.1 Standard light curve 



For convenience, we put the lens at the origin, and let the source move across the line of sight 
along the x-axis (see Fig. The impact parameter in units of the Einstein radius is labelled as uq. 
For convenience, we define the Einstein radius crossing time (or 'timescale') as 



'E — — — , t/E — — , /Irel = — 

Vt A^rel Dd Dd 



(4.1) 



where Vt is the transverse velocity and /^rei is the relative lens-source proper motion. Substituting 



the expression for the Einstein radius into eq. (3. 1 1 ), we find that 



Vt 



.E-19dayJ4xg(l-g) (^) ^ ^ V200kms-i 



-1 



(4.2) 



If the closest approach is achieved a time t = to, then the dimensionless coordinates are Xj = (f 
to)/tE and = Mo, and the magnification as a function of time is given by 



nit) 



t-ta 



r,{t)^r,^{t)+4 - V " V ?E 
Two light curve examples are shown in the right panel of Fig. |l| for mq = 0.1 and 0.3. 



(4.3) 



* Source 



Lens 



Figure 4: Illustration for the lens position and source trajectory. The dimensionless impact parameter is uq, 
{xs,ys) are the dimensionless source position along the trajectory, and rs is the distance between the lens and 
source. 



To model an observed light curve, three parameters are present in eq. ( [4.3[ ): to, tE, mq- In 
practise, we need two additional parameters, mo, the baseline magnitude, and /s, a blending pa- 
rameter, /s characterises the fraction of light contributed by the lensed source; in crowded stellar 
fields, each observed 'star' may be a composite of the lensed star, other unrelated stars within the 
seeing disk and the lens if it is luminous ([^, 32 1). Blending will lower the observed magnification 



and in general /s depends on the wavelength, and so each filter requires a separate /s parameter. 



Unfortunately, we can see from eq. (4.3) that there is only one physical parameter (fg) in the model 



that relates to the lens properties. Ie depends on the lens mass, distances to the lens and source, 
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and the transverse velocity /Xrei- Therefore from an observed Ught curve well fitted by the standard 
model, one cannot infer the lens distance and mass uniquely; this is the so-called microlensing 
degeneracy. However, given a lens mass function and some kinematic model of the MiUcy Way, we 
can infer the lens mass statistically. 

4.2 Non-standard light curves 

The standard model assumes the lensed source is point-like, both the lens and source are single 
and all the motions are linear. The majority (~ 90%) of microlensing events are well described by 
this simple model. However, about 10% of the light curves are non-standard (exotic), due to the 
breakdown of one (or more) of the assumptions. We briefly list these possibilities below (see the 
talk by Dominik for more details.) These non-standard microlensing events allow us to derive extra 
constraints, and partially lift the microlensing degeneracy. Because of this, they play a role far 
greater than their numbers suggest. 



(1) The lens may be in a binary or even a multiple system ([|35[]). The light curves for a binary 
or multiple lensing system can be much more diverse (see i3). They offer an exciting way 
to discover extrasolar planets (||3|, |7[ |[ ||]). 

(2) The source is in a binary. In this case, the light curve will be a simple, linear superposition 
of the two sources (when the rotation can be neglected, see [p^]). 

(3) The finite size of the lensed star cannot be neglected. This occurs when the impact parameter 
Mo is comparable to the stellar radius normalised to the Einstein radius, mq ^ r^^/r-^. In this 



case, the light curve is significantly modified by the finite source size effect ([]56|, [23|]). The 
finite source size effect is most important for high magnification events. 

(4) The standard light curve assumes all the motions are linear. However, the source and/or the 
lens may be in a binary, furthermore, the Earth rotates around the Sun. All these motions 
induce accelerations. The effect due to the Earth motion around the Sun is usually called 
"parallax" (e.g. [^, ^]), while that due to binary motion in the source plane is called 
"xallarap" ("parallax" spelt backwards, [Q, |I|]). Parallax or "xallarap" events usually have 
long timescales. For a typical microlensing event with timescale ~ 20day, the parallax 
effect due to the Earth rotation around the Sun is often undetectable (unless the photometric 
accuracy of the light curve is very high). 

(5) Microlensing can "repeat", in particular if the lens is a wide binary ([|l5|]) or the source is a 
wide binary. In such cases, microlensing may manifest as two well-separated peaks, i.e., as 
a "repeating" event. A few percent of microlensing events are predicted to repeat, consistent 
with the observations (||50|]). 

Notice that several violations may occur at the same time, which in some cases allow the microlens- 
ing degeneracy to be completely removed (e.g. 1^, 18, 21 1). 
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4.3 A'^-point lens gravitational microlensing 

It is straightforward to derive the (dimensionless) lens equation for A/^-point lenses. We can 



first cast eq. (3.8) in vector form and then rearrange 

^ ^ 1 



r. (4.4) 



The above expression implicitly assumes that the lens is at the origin, and all the lengths have been 
normalised to the Einstein radius corresponding to its mass (or equivalently, the lens mass has been 
assumed to be unity). 

Let us consider the general case where we have A'^-point lenses, at = {xk,yk) with mass M^, 
k= I, - ■ - jN. We normalise all the lengths with the Einstein radius corresponding to the total mass, 
M = Y!k=i ^k- The generalised lens equation then reads 

ti \{r-rk)V M 



where Lf^i = 1. If there is only one lens (mi = 1) and the lens is at the origin, then eq. ( |4-.5D 
reverts to the single lens equation (p^. 

Two-dimensional vectors and complex numbers are closely related, Witt (1990) first demon- 
strated that the above equation can be cast into a complex form by direct substitutions of the vectors 
by complex numbers: 

V' ^ ~ ^ m , . 

i^l {Z-Zk){z-Zk) k^lZ-Zk 

where z = x + yi,Zk=Xk +yk U and Zs =x^+y^ i (where / is the imaginary unit). 



We can take the conjugate of eq. (4.6) and obtain an expression for z- Substituting this back 
into eq. ( |4.6| ), we obtain a complex polynomial of degree N'^ + \. This immediately shows that 
even a binary lens equation cannot be solved analytically since it is a fifth-order polynomial ^. 

The magnification is related to the determinant of the Jacobian of the mapping from the source 



plane to the lens plane: (jCg, js) {x,y). In the complex form, this is ([]55|]; see also the solution to 
Problem 8): 

^^1^ j^ d{x,,y,) ^^ dz^'d^^ ^4^^ 
^(■^,3') dz dz 

Notice that the Jacobian can be equal to zero implying a (point) source will be infinitely magnified. 
The image positions satisfying 7 = form continuous critical curves, which are mapped into caus- 
tics in the source plane. Of course, stars are not point-like, they have finite sizes. The finite source 
size of a star smoothes out the singularity. As a result, the magnification remains finite. 



For A/^-point lenses, from the complex lens equation (^.6), we have 

dzs ^ m , , 



dz H{z-Zkf' 



nik 



k^\ {z-Zkf 



(4.8) 



^In classical mechanics, the two-body problem can be solved analytically, but not the three-body problem. 
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Figure 5: Left panel: Caustics (red curves) and critical curves (green curves) for a binary lens. The lenses 
(indicated by two '+' signs) are located at zi = (—0.665,0) and Z2 — (0.035,0) with mass nii — 0.95 and 
and OT2 = 1 — nil = 0.05 respectively. The black line shows the trajectory for three source sizes, rs/r^ — 
0,0.05,0.3, indicated by the cyan and blue circles and a dot (for a point source). The trajectory starts at 
(—2,-1) with a slope of 0.7. Right panel: Corresponding light curves for the three source sizes along 
the trajectory in the left panel. The time is normalised to the Einstein radius crossing time, fg, and t — 
corresponds to the starting position. Notice that as the source size increases, the lensing magnification 
amplitude decreases. 



It follows that the critical curves are given by 



1. (4.9) 



The sum in the above equation must be on a unit circle, and the solution can be cast in a parametric 
form 

i=l {z-ZkY 

where < < 2?! is a parameter. The above equation is a complex polynomial of degree of 2N 
with respect to z- For each there are at most 2N distinct solutions. As we vary <I> continuously, 
the solutions trace out at most 2N continuous critical curves (critical curves of different solutions 
may join with each other smoothly). In practise, we can solve the equation for one value, and 
then use the Newton-Raphson method to find the solutions for other values of <I>. 

For a single point lens, if we take zi = 0, and m\ = 1, we find that the critical curve is the 
Einstein ring (|z| = 1), which is mapped into a degenerate caustic point at the origin (zs = 0). 
However, for binary or multiple lenses, the critical curves and caustics are much more complex. 
The left panel in Fig. || illustrates the critical curves and caustics for a binary lens with mi = 0.95 
and m2 = 0.05 and separation of 0.7 (in units of the Einstein radius for the total mass). In this case, 
there are three separate, continuous critical curves which are mapped into three caustics. 

For a point source, the complex polynomial can be easily solved numerically (e.g. using the 
z roots routine in [ff|]). However, for a source with finite size, the existence of singularities 
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makes the integration time-consuming (see section The right panel in Fig. ^ shows the light 
curves for three source sizes along the trajectory indicated by the green straight line. As the source 
size increases, the lensing magnification amplitude decreases, and the number of peaks differ for 
different source sizes. 

5. Optical depth and event rates 

So far we have derived the lens equation and light curve for microlensing by a single star. 
In reality, hundreds of millions of stars are monitored, and 800 unique microlensing events 
are discovered each year. Clearly we need some statistical quantities to describe microlensing 
experiments. For this, we need two key concepts: optical depth and event rate. 

5.1 Optical depth 

The optical depth (lensing probability) is the probability that a given source falls into the 
Einstein radius of any lensing star along the line of sight. Thus the optical depth can be expressed 
as ^ 

r = J^ \{Da) {KrE^) dDa, (5.1) 

which is an integral of the product of the number density of lenses, the lensing cross-section (= 
TTrn^) and the differential path (dDd). 

Alternatively, the optical depth can be viewed as the fraction of sky covered by the angular 
areas of all the lenses, which yields another expression 

T = ^ £^ [n{Da)4KDd^dDa] {ndE^) , (5.2) 

where the term in the [ ] gives the numbers of lenses in a spherical shell with radius to + dDd, 
kOe' is the angular area covered by a single lens, and the term in the denominator is the total solid 
angle over all the sky {An). 

If all the lenses have the same mass M, then n{Di) = p(Dd) /M, nr-^' °^ M, and the lens mass 
drops out in n{Dd)7lrE^. Therefore the optical depth depends on the total mass density along the 
hne of sight, but not on the mass function. 

Let us consider a simple model where the density is constant along the hne of sight, p (Dd) = 
Pq. Integrating along the line of sight one finds 

_ 2nGpo ^2_ 1 Gpo47iD,y3 _ 1 GM{<D,) _V^ 

where M{<D^) is the total mass enclosed within the sphere of radius and the circular velocity 
is given by = GM{< D,) /D,. 

For the Milky Way, V 220kms^^ t ps 2.6 x 10^^. The low optical depth means millions 
of stars have to be monitored to have a realistic yield of microlensing events, and thus one needs 
to observe dense stellar fields, which in turn means accurate crowded field photometry is essential 
(see the talk and workshop material by P. Wozniak on difference image analysis^.). 

**pos . sissa. it//archive/conferences/054/003/GMC8_003 .pdf ; 
www . jodrellbank . manchester .ac.uk/~m2c/ dia_exercises/ 
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5.2 Event rate 

The optical depth indicates the probability of a given star that is within the Einstein radii of 
the lenses at any given instant. As such, the optical depth is a static concept. We are obviously 
interested in knowing the event rate (a dynamic concept), i.e., the number of (new) microlensing 
events per unit time for a given number of monitored stars, A^^. 

To calculate the event rate, it is easier to imagine the lenses are moving in a static stellar source 
background. For simplicity, let us assume all the lenses move with the same velocity of Vt. The 
new area swept out by each lens in the time interval dt is equal to the product of the diameter of 
the Einstein ring and the distance travelled Vtdt, dA = 2rE x vidt = Ir-^'dt /t-E- The probability of a 
source becoming a new microlensing event is given by 

di = j n{Dd)dAdDd = J n{Di) i-j- \ dtdD^ (5.4) 
The total number of new events is N^^d T, and thus the event rate is given by 

r= ^(M /°-„(o,) = f'-l. (5^5) 

dt Jo yntE J 71 J tE 

If, for simplicity, we assume all the Einstein radius crossing times are identical, then we have 



K tE 

We make several remarks about the event rate: 



(5.6) 



(1) If we take tE = 19day (roughly equal to the median of the observed timescales), then we 
have 

r~— - = 1200yr-^— ^— ^ — i— (5.7) 
K tE 10^ 10-6 \^l9day J ^ ' 

For OGLE-III, about 2 x 10^ stars are monitored (see Udalski's contribution), so the total 
number of events lenses we expect per year is F ~ 2400 if T ~ 10^^, which is a factor of four 
of the observed rate (indicating the detection efficiency may be of the order of 30%). 

(2) While the optical depth does not depend on the mass function, the event rate does because of 
fE(oc M^/^) in the denominator of eq. (^I^. The timescale distribution can be used to probe 
the kinematics and mass function of lenses in the Milky Way. 

(3) The lenses and sources have velocity distributions, one must account for them when realistic 
event rates are needed. Furthermore, the source distance is unknown, and so in general we 
need to average over the source distance (for example calculations, see [ ^ , |1|]). 



6. Summary 

In this introduction, we derived the lens equation, and obtained the image positions and magni- 
fications for a point lens. We also discussed the statistical measures for microlensing experiments, 
and estimated the order of magnitudes for various quantities. An interested reader should now be 
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armed with a basic knowledge of microlensing and be prepared to read the review articles and start 
to do research on gravitational microlensing (or even try to solve the problem set below). 

Since the discovery of first microlensing events in the early 1990's, enormous achievements 
have been made in the field. However, challenges and opportunities remain. 

(1) Undoubtedly the highlight of gravitational microlensing in the last few years has been the 



discovery of extrasolar planets ([|12|, p3|, a, 24, 11]). Microlensing has much to offer in this 



area since it probes a different part of the parameter space, and provides an important test 
of the core accretion theory of planet formation. Several White Papers ([25, 10, [l7|, ^) set 



out strategies with ambitious milestones in the next fifteen years, from improvement of the 
current survey plus foUowup mode of discovery (with an automated algorithm to identify 
the "anomalies" in real-time) in the near term, to a wide-field network from the ground in 
the next 5-10 years, and eventually a telescope in space in the next 10-15 years. Combined 
with the stellar transit mission Kepler (to be launched in 2009), microlensing will be able to 
provide the complete census of Earth-mass (and lower) planets at virtually all the separations. 

Technically, it is still challenging to calculate the light curves for sources with finite size 
since we need to integrate over the singularities of caustics. This is particularly important 
for the discovery of extrasolar planets when a source transits the small caustics induced by 
the planet(s). The problem becomes even worse with the discovery of multiple planets ([^]) 
due to the higher complexity of the lens equation and the increased number of parameters: 
how do we search the high dimensional parameter space efficiently? 

Are there hidden multiple planetary light curves in the database that are not yet identified 
due to their complex shapes? 

(2) Microlensing surveys over the last fifteen years have accumulated tens of TB of data. This 
tremendous database has not been exploited to its fullest potential. 

For example, the surveys yielded many high-quality colour-magnitude diagrams of stellar 
populations, proper motions of millions of stars, and in the future the optical depth maps. 
All these can be used to provide important and independent probes of the structure of the 
Milky Way. 



Despite promising earlier attempts (e.g. [37]; [M, 45]; [31, pq]), microlensing has under 



delivered in this area. For example, while we have discovered several thousands of mi- 
crolensing events over the last 15 years, only a small fraction has been used for statistical 
analyses of optical depths. We need to remedy the situation urgently. 

(3) High-magnification events are great targets-of-opportunity for high signal-to-noise ratio spec- 
troscopic observations to study stellar atmospheres for bulge stars. Attempts so far already 



yielded interesting results (e.g. [63, 52, 13]). We need to explore this more systematically. 



(4) For mathematically-gifted students (or mathematicians), gravitational microlensing provides 
an interesting problem. While the binary lens equation is no longer analytical, there is, nev- 
ertheless, an analytical relation on the minimum magnification for five-image configurations 
([^, 47]). There is also a degeneracy found by Dominik (1999) between close and wide 
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separation binaries, which was later explored in much greater detail by An (2005). Are there 
any other symmetries, perhaps even for multiple lens systems? 



The number of critical curves for A/^-point lenses cannot exceed 2N (see section 4.3). The 
upper bound of the number of images for A/^-point lenses also has linear dependence on 
(see Problem 1). These two are related mathematically (see [48, 49, ^]). Is there a more 
geometric (topological) way of understanding this? 



I thank Eamonn Kerins, Szymon Kozlowski, Nick Rattenbury and particularly Matthew Penny 
for a careful reading of the manuscript. 
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Problems 

1. a) There are A^-point lenses in a lens plane. For a source very far away from all the lenses, 
how many images are there? What are their parities? 

b) What is the achievable maximum number of images for A^-point lenses and an arbitrary 
source position? 

2. A uniform source star is perfectly aligned with an observer and a point lens. Its physical 
radius normahsed to the Einstein radius is p*. 

a) What is the resulting image configuration? 

b) What is the magnification for the finite source? 

c) Estimate the maximum magnification that can be achieved for a source star in the Galac- 
tic bulge. 

d) Derive the expression for magnification when the source is not perfectly ahgned with 
the lens. 

3. Show that the total magnification for a point lens is always larger than one. How can this be 
reconciled with energy conservation? (see Jaroszyhski & Paczynski 1996, AcA, 46, 361). 

4. The density distribution in the plane of the Galactic disk can be modelled as an exponential 

p{R)=PoQxip{-{R-Ro)/Rd), 

where po is the density in the solar neighbourhood, Rq is the distance from the Sun to the 
Galactic centre, R^ is the disk scale length, and R is the distance from the Galactic centre to 
the lens. 

a) Find the optical depth T for a source at the Galactic centre {R = 0). 

b) If Po = 0. \Mq pc~^, /?o = 8 kpc, Rd = ^ kpc, what is the value of T? 
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5. Consider a simple model: all lensing objects have the same mass M, the same three-dimensional 
velocity V, and their velocity vector directions have an isotropic distribution. The source 
located at the distance is stationary, and the number density of lensing objects is uniform 
between the observer and the source. Derive the timescale probability distribution. 

Now assume the lenses follow a Maxwellian distribution with a one-dimensional velocity 
dispersion a. Derive the timescale probability distribution. Show that it follows a power-law 
behaviour for both very short and very long timescales. 

6. A distance source is lensed by a point deflector with mass M. The light signals emitted by 
the source will be received by an observer at different times for the two images due to the 
difference in the trajectory and gravitational potential experienced. The time delay is of the 
order of rjch/c, where r^ch is the Schwarzschild radius. Is this observable for a M = IMq 
lens? 

7. A background star stationary at the origin is microlensed by a lens moving from — oo to o°. 
Show that the centre of light of the two images traces out an ellipse. 

8. Show that 



a) In the complex notation, the Jacobian is given by eq. (|4.8|). 

b) For any positive-parity image produced by A^-point lenses, the magnification is always 
larger than or equal to one. 

c) Find the number and positions of images where their magnifications are identical to 
unity for a binary lens. 

d) What is the maximum number images with unity magnification for A/^-point lenses? 
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Solutions 



a) There are A'^ + 1 images, one bright image close to the source, and A'^ fainter images 
with one image close to each lens. The image close to the source has positive parity, all 
the others have negative parities. 

b) 5(A^ - 1) for N>2. This is a difficult problem. See Rhie S. H. 2003, arXiv:astro- 
ph/0305166; Khavinson, D. and Neumann G. 2006, Proceedings Of the American 
Mathematical Society, 134, 1077 (arXiv:math/0401 1088). 



2. a) The image configuration is an annulus. The outer radius is given by the r2 = ( + 

p*)/2 while the inner radius is given by n = (Vp* +4 — P*)/2. 

b) The area covered by the image is A = Tz{r^ ~ 'i) = ^p*\/p| + 4. Thus the magnifica- 
tion is n =A/{kpI) = vW+^/P*- In particular, when p^^O, ^ 2/p*. 

c) For a solar mass lens in the Galactic centre (D^ = 8 AU), the maximum Einstein radius 
is achieved when = Df^/l, = 6.0 x 10^^ cm. The faintest and smallest stars we 
can see in the bulge have radii similar to a solar- type star, f» /?o = 7 x lO^^'cm, p^ = 
r^/rE = 1.2 X 10^1 The maximum magnification is of the order of ml/pi^m 1700. 

d) See Witt & Mao 1994, ApJ, 430, 505 

3. Clearly > 1. See Jaroszynski & Paczynski 1996, AcA, 46, 361 for discussions about the 
energy conservation. 

4. The optical depth is given by 

Jo "(^d) (^^ ) = -2 ^^^d- (6.1) 

For microlensing in the Galactic plane, Djs = R,Dfi = Ro^Dd =Ro—R, and thus we have 

po e"^/*" dDd (6.2) 



DJR, 47tG Dd{Ro- Dd) 
c2 Ro 



Performing the above integral, we find that 



z = ^^^^y-'[2 + y + ey{-2 + y)], y = Ro/Rd (6.3) 



For the given numbers, we find4;rGpo/?§/c2 = 3.86 x 10"^ andy = /?o/^rf = 2.67, we have 



2.9 xlO"^ (6.4) 



5. For step by step derivations, see Mao & Paczynski (1996, ApJ, 473, 57). Notice that the 
observed event timescale distribution does not follow the power-laws due to detection effi- 
ciency. 

6. The time delay is of the order of tens of micro-seconds for a solar mass lens, and is very 
difficult to observe. 



18 



Introduction to Gravitational Microlensing 



Shude Mao 



7. See Hog et al. (1995, A&A, 294, 287), Miyamoto & Yoshi (1995, AJ, 110, 1427), Walker 
(1995, ApJ, 453, 37) 



8. a) Using complex notation, we have x = ^{z + z) ,y = j]{z — ■ Thus 



dZs _ dzs_d^ _^ dz^dy^ _ 1_ / dzs .dzs 
dz dx dz dy dz 2 \ dx dy 



dx^ ^ .dys _ .dx^ _^ dy^ 
dx dx dy dy 



(6.5) 



Similarly 

dzs _ dzs dx ^dzs dy _ \ / dzs _^ .dZs 
dz dx dz dy dz. 2 \ dx dy 



1 / dxs . dy. 



-^a , 1 - / — » , .-J!, . .dx^ dy^ . 



However, from the lens equation ( H^ , dzs/dz = 1, and thus comparing the real and 

dx. 



imaginary parts in eq. (6^), we have 



1 / ^ dy, 

2 \dx dy 



1, 



dys 

dy dx 



Substituting the second expression into eq. (|6^), we find 

dzs 
dz 



2 \ dx dy J dx 



Combined with 



dZf. 



I f dxfi ^ dy, 
dx dy 



(6.7) 



(6.8) 



(6.9) 



we find 

dx^ _ 
dx 



1 + 



1 / 5zs , dzA dys 



■ + 



dz dz j dy 



1 



I ( dZf. , dZfi 



dz ~'~ dz 



dys 
dx 



1 / dz. 



2i \ dz 



dz, 

dz 



Substituting the above equations into the Jacobian 

j_ dx, dy, dy, dx, 
dx dy dx dy ' 

we recover the required expression, 
b) Since 

7=1 



(6.10) 



(6.11) 



(6.12) 



For a positive parity image, we must have 1 > 7 > 0, it follows that jx = J > 1 . The 
magnification is unity when the sum is equal to zero. 

c) Without losing generality, we put the lenses on the ;c-axis, at r\ = {zo,0) and r2 = 
(— Z0)0), with masses m\ and m2 {m\ > 022), and we have m\+m2 = 1, Zo > 0- From 
part b), the condition for unity magnification is 

nil fni 



+ ■ 



(z + zo)^ (z-zo)2 



0. 



(6.13) 



19 



Introduction to Gravitational Microlensing 



Shude Mao 



This simplifies into a quadratic equation 

^ + 2{mi-m2)zoz + zi = 0. (6.14) 
Clearly we have two solutions 

z+- = ZQ (^m2 - mi) ± isj 1 - (mi - mi)^^ . (6.15) 

These two are, as expected, symmetric with respect to the;c-axis. 

d) 2{N-l). 
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